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ABSTRACT 


At  low  temperatures  the  vibrational  relaxation  time,  t  ,  is 
much  shorter  than  the  dissociation  time,  t^.  The  -  Ar  results  of 
Camac  (J.  Chem.  Phys.  34  448  (1961)  )  and  Camac  and  Vaughan 
(Ibid,  p.460)  yield  a  T(j/Tv  =  60  at  5000°K  and,  upon  extrapolation, 
a  T(j/Tv  =  1.  4  at  18000°.  According  to  these  extrapolations,  dissociation 
at  high  temperatures  would  proceed  significantly  before  vibrational  equili¬ 
bration  would  occur.  The  purpose  of  this  investigation  was  to  determine 
how  the  dissociation  rate  will  be  affected  by  a  lack  of  vibrational  equilibrium. 

Studies  of  the  dissociation  rate  of  dilute  -  Ar  mixtures  were 
made  in  a  24"  diameter  shock  tube  from  5000  -  18000°K.  The  con¬ 
centration  was  monitored  by  its  absorption  of  1470A  radiation.  An 

Arrhenius  plot  of  the  data  yielded  a  straight  line  from  5000  -  11000°K, 

14 

the  rate  constant  being  given  by  k^  =  2.  9  (+  12%)  x  10  exp  (-D/RT) 
cc/mole-sec.  Above  11000°  the  data  deviate  from  the  line  given  by  this 
equation  -  at  18000°  k^  being  .  45  times  the  calculated  value. 

An  incubation  time,  At,  was  observed  during  which  dissociation 
does  not  proceed  to  a  significant  extent.  The  ratio  of  this  incubation 
time  to  the  vibrational  relaxation  time  (obtained  by  extrapolating  Camac' s 
low  temperature  results)  when  plotted  against  translational  temperature, 
displays  a  slight  negative  temperature  dependence.  At  18000°  At/rv  =  0.  4, 
at  8000°  At ;/t  =  1,  and  at  5500°  we  estimate  that  At/r^  w  2. 
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I.  Introduction 

The  rates  of  dissociation  and  vibrational  relaxation  of  oxygen  for 
the  catalytic  species  Ar,  O2  and  O  are  now  reasonably  well  determined  up 
to  T  =  7500°K.  We  quote  here  for  the  rate  of  dissociation  of  oxygen  by 
argon  the  results  of  Camac  and  Vaughan* 

^=3.6x10*®  T~*  exp  (-D/RT)  cc/mole-sec  ,  (1) 

where  D  =  118000  cal/mole.  A  dissociation  relaxation  time,  r  is 
defined  by 

*Ad*kd  (nAr/N0  )  •  <2> 

where  n  .  is  the  number  of  argon  atoms/cc  and  Nn  is  Avogadro's  number. 
Ar  ”  2 

For  the  vibrational  relaxation  time,  r  ,  of  by  Ar,  Camac 

gives  r  -1 

1/tv=1.  2x10'7  T1/6  I  l-exp(-2228/T)  I  exp(-1.04xl07/Tr  3nArsec-1.  (3) 

v  $ 

These  equations  are  plotted  in  Fig.  1  setting  n  =1.28  x  1017. 

At  low  temperatures  the  vibrational  relaxation  process  is  much 
faster  than  the  dissociation  process.  Hence,  the  rate  constant  which  is 
determined  for  the  latter  process  is  for  a  gas  essentially  in  vibrational 
equilibrium.  As  the  temperature  is  increased  the  rates  of  both  processes 
increase  but  the  dissociation  rate  increase-s  faster;  according  to  Fig.  1, 
at  T  =  5000°K,  r  Jt  y  =  60;  at  T  =  8000°,  t  ^/t  v  =5.4  and  at  T  =  18000°, 
Td/Tv  =  1-  4.  These  ratios  are  independent  of  density.  Hence,  at  tempera¬ 
tures  above  about  8000°,  on  the  basis  of  the  above  extrapolations,  one 
would  expect  that  dissociation  would  proceed  to  a  significant  extent  before 

^  This  particle  density  corresponds  to  T  =  300°K  and  P  =  4mm  (or  an 
initial  shock  tube  pressure  of  P.  =  1mm  and  a  density  ratio  across  the 
shock  of  Po/Pj  =  4).  This  density  was  chosen  as  it  corresponds  to  the 
♦  experimental  conditions  used  in  this  study. 
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Fig.  1  The  dissociation  and  vibrational  relaxation  times  in  the  particle 
coordinate  system  at  a  density  of  n.  =  1.  28  x  10^  particles/cc 
plotted  as  a  function  of  temperature:  The  T&  is  extrapolated  from 
experimental  results  above  7500°,  tv  above  7000°.  Also  indicated 
in  the  figure  is  the  density  independent  ratio  t(j/tv  at  several 
temperatures. 


vibrational  equilibrium  obtains.  The  purpose  of  this  experimental  study 
«  was  to  determine  how  the  dissociation  rate  will  be  affected  by  this  situation, 

i.  e.  ,  to  what  extent  the  dissociation  and  vibrational  processes  couple. 

Camac  did  his  work  in  a  1.  5  inch  diameter  shock  tube  where  it 
was  not  possible  to  obtain  adequate  resolution  at  temperatures  at  which 
coupling  of  vibration  and  dissociation  might  become  important.  However, 
by  using  a  21.  5%  -  78.  5%  Ar  mixture  at  an  initial  pressure  of  0.  93mm 

he  did  observe  a  coupling  effect  at  temperatures  around  9000  .  It  should 
be  noted  that  with  such  a  rich  mixture,  there  are  large  temperature 
changes  behind  the  shock  front  as  the  relaxes;  this  tends  to  complicate 
the  data  analysis.  Camac  states  that  his  data  show  that  dissociation  "is  in 
progress  during  the  vibrational  relaxation"  although  it  is  slower  "by  more 
than  a  factor  of  2". 

In  the  present  work,  the  dissociation  of  dilute  O^-Ar  mixtures  has 
been  measured  using  a  24"  diameter  shock  tube  available  at  this  labora¬ 
tory.  3  This  facility  allows  the  production  of  shock  waves  at  the  low  ini¬ 
tial  densities  necessary  to  slow  down  the  relaxation  processes  sufficiently 
to  be  well  resolved.  Furthermore,  the  long  optical  path  allows  a  low 

* 

partial  pressure  of  O 2,  hence,  essentially  isothermal  and  isobaric  con¬ 
ditions  behind  the  shock.  The  temperature  range  covered  in  this  investi¬ 
gation  was  5000  -  18000°K.  The  C>2  concentration  was  monitored  as  a 
function  of  time  behind  the  incident  shock  wave  by  its  absorption  of  1470A 

radiation  -  similar  to  the  work  of  Camac. 

2. 

Camac  has  measured  the  absorption  coefficient  of  oxygen  at 
1470A  as  a  function  of  its  vibrational  temperature  and  showed  that  this 
absorption  coefficient  is  independent  of  the  rotational  and  translational 
temperature  of  the  gas.  At  room  temperature  the  absorption  coefficient 
per  cm  of  path  length  at  a  density  corresponding  to  300°  and  one  atmos¬ 
phere  pressure  is  315.  He  found  that  this  coefficient  decreases  mono- 
tonically  as  the  vibrational  temperature  is  increased,  and  asymptoti¬ 
cally  approaches  the  value  of  135  at  about  7500°  and  remains  constant 
thereafter.  Although  his  measurements  of  this  absorption  coefficient 
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were  carried  just  a  little  past  8000°,  theoretical  arguments  justify  assum¬ 
ing  that  it  remains  constant  as  one  continues  to  increase  the  temperature. 
These  results  are  plotted  in  Fig.  2. 

At  the  high  end  of  the  temperature  range  covered  in  the  present 
work  the  equilibrium  ionization  of  the  argon  diluent  is  very  significant  and 
the  energy  invested  in  ionization  would  cause  the  temperature  to  drop 
many  thousands  of  degrees.  In  order  to  ascertain  whether  or  not  ioniza¬ 
tion  did  occur  during  the  vibrational  and  chemical  relaxation  of  the  small 
amount  of  in  the  gas,  the  electron  density  was  monitored  using  a 
microwave  technique.  ^ 

II.  The  Apparatus 

-  2 

The  shock  tube  used  in  this  work  was  designed  and  built  by  Lin 

for  his  investigations  of  the  ionization  rate  in  shock -heated  air,  although 
for  the  present  work  a  new  test  section  was  built.  The  tube  is  of  stainless 
steel  construction  having  a  length  of  60'  and  a  diameter  of  2'.  The  driver 
section  is  6'  long  and  the  diaphragm  is  ruptured  by  igniting  mix¬ 

tures  at  high  pressure.  Shock  speeds  were  determined  using  painted 
platinum  strip  heat  transfer  gages. 

Because  there  was  some  concern  over  the  possibility  of  the  shock 
wave  being  tilted  with  respect  to  the  shock  tube  walls,  four  heat  transfer 
gages  were  mounted  around  the  circumference  of  the  shock  tube  90°  apart 
just  in  front  of  the  optical  station.  The  signals  from  two  of  these  gages 
were  delayed  three  microseconds  and  put  into  the  -B  channels  of  the 
differential  amplifiers  of  a  dual  beam  scope.  The  two  undelayed  signals 
were  put  into  the  +A  channels.  In  this  way  we  could  correlate  the  time  of 
arrival  of  the  shock  wave  at  the  four  heat  transfer  gages.  Under  no  condi¬ 
tions  in  these  studies  was  the  shock  front  tilt  found  to  be  significant  (i.  e.  , 
it  was  less  than  the  optical  slit  width). 

At  low  densities  a  shock  front  in  a  shock  tube  is  curved  to  a 
3 

significant  extent.  Lin  has  measured  the  shock  front  curvature  as  a 
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function  of  initial  pressure  and  shock  speed  for  air  shocks  in  this  24"  dia¬ 
meter  tube.  He  found  the  shock  front  curvature  (defined  as  the  distance 
along  the  center  line  from  the  shock  front  to  the  plane  formed  by  the 
intersection  of  the  shock  front  and  the  shock  tube  walls)  to  be  independent 
of  the  shock  speed  and  to  vary  inversely  as  the  square  root  of  the  initial 
pressure.  His  data  for  curvature,  6,  as  a  function  of  initial  pressure  P^, 
can  be  fit  by  the  following  equation: 

6  (mm)  =  1.  75  (mm  of  Hg)  .  (4) 

At  the  lowest  initial  pressures  used  in  these  experiments  (1/4 
mm),  the  shock  front  curvature  would  amount  to  3.  5  mm.  The  optical 
slits  used  for  all  the  experiments  were  1  mm  wide,  hence,  a  shock  front 
curvature  of  3.  5  mm  represents  a  serious  loss  in  resolution.  To  avoid 
this  difficulty  the  optical  path  was  reduced  to  12"  by  inserting  1"  diameter, 

6"  long  pipes  into  opposite  sides  of  the  shock  tube.  (See  Fig.  3.  )  Each  of 
these  pipes  was  terminated  inside  the  shock  tube  by  a  flat  plate  2.  5"  high 
and  1.  5"  wide,  the  edges  of  which  were  tapered  at  an  angle  of  9°  44'  to 
form  knife  edges.  *  The  center  of  each  knife-edged  plate  contained  a  flush 
calcium  fluoride  window  1/2"  in  diameter  and  3/8"  thick;  behind  each 
window,  i.e.  ,  inside  the  pipe  was  located  the  1/2"  high,  1mm  wide  slit. 
Assuming  the  shock  front  surface  to  be  spherical,  the  amount  of  curvature 
between  two  parallel  planes  placed  within  the  shock  tube  decreases  as 
the  square  of  the  ratio  of  the  distance  between  the  two  planes  to  the  shock 
tube  diameter.  Hence,  for  the  l/4  mm  runs  the  shock  front  curvature  was 
effectively  reduced  to  .87  mm.  Furthermore,  this  technique  reduces  the 
problem  of  0-atom  recombination  in  the  cold,  dense  boundary  layer  since  on 


*  The  width  of  the  plate  was  arrived  at  by  making  the  leading  edge  project 
twice  the  standoff  distance  of  the  bow  shock  formed  in  front  of  the  pipe 
and  the  height  was  such  that  the  Mach  wave  formed  at  the  upper  and  lower 
corners  of  the  leading  edge  would  not  cross  the  window  located  at  the  center 
of  the  plate. 
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Fig.  3  Schematic  diagram  of  the  shock  tube  test  section.  Both  the  U.  V. 
optical  system  and  the  microwave  system  are  shown. 


these  plates  the  boundary  layer  grows  only  to  a  steady  state  thickness 
instead  of  growing  continuously  as  on  the  shock  tube  wall. 

The  light  source  was  located  externally  to  the  shock  tube  at  the 
entrance  to  one  of  these  pipes  (see  Fig.  3).  The  light  source  was  a 
Lyman  lamp,  which  produced  an  oscillating  output  whose  second  cycle 
maximum  was  flat  for  about  20  microseconds  (see  Ref.  2  for  a  description 
of  this  light  source ).  This  lamp  was  triggered  from  one  of  the  upstream 
heat  transfer  gages  through  a  time-delay  generator  which  was  pre-set  to 
synchronize  the  second  cycle  light  maximum  with  that  of  the  shock  wave 
arrival  at  the  optical  station.  The  monochromator,  consisting  of 
two  calcium  fluoride  prisms  and  necessary  light  baffles,  was  located 
outside  the  shock  tube  at  the  entrance  to  the  other  pipe.  A  photomultiplier 
coated  with  sodium  salicylate  phosphor  was  mounted  within  the  mono¬ 
chromator  just  behind  the  exit  slit.  In  use,  the  monochromator  was 
evacuated.  The  output  of  the  photomultiplier  was  monitored  on  an  oscillo¬ 
scope.  This  oscilloscope  was  triggered  by  a  heat  transfer  gage  located 
about  one  inch  in  front  of  the  optical  station. 

The  electron  density  monitor  was  mounted  perpendicular  to  the 
U.  V.  system.  In  this  case  the  wave  guide  itself  served  as  the  six  inch 
pipe  extending  into  the  shock  tube.  It  was  also  terminated  by  a  knife -edged 
plate  similar  to  that  of  the  U.  V.  system,  the  open  end  of  the  wave  guide  being 
covered  by  an  oblong  quartz  window  (height  =  .  311",  width  =  .  622").  The 
microwave  frequency  used  was  f  =  1.  5  x  10^  cps  (Ku  band).  Before  each 
run,  the  maximum  deflection  of  the  oscilloscope  trace  from  the  reflected 
microwave  signal  was  set  by  shorting  out  the  wave  guide  (corresponding 
to  complete  reflection).  This  microwave  reflection  signal  was  displayed 
on  a  dual  beam  scope,  the  other  beam  of  which  displayed  the  U.  V.  signal. 


III.  Data  Reduction 

Three  different  composition  C^-Ar  mixtures  were  used,  they 
were:  4%  02  (Pj  ~l/4  mm),  2%  Oz  (P1  ~  1  and  1/2  mm)  and  1/2% 

(P^  ~4  and  10  mm).  The  richer  mixtures  were  employed  at  higher 
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temperatures  where  lower  densities  were  necessary  to  allow  resolution. 

At  lower  temperatures,  where  the  rate  processes  are  slower,  higher 
densities  could  be  used  in  conjunction  with  smaller  O2  partial  pressures. 
This  is  a  desirable  situation  since  the  relaxation  processes  are  most 
temperature  dependent  at  low  temperatures. 

For  each  of  the  three  different  mixtures  studied,  shock  tube 
performance  data  were  computed  and  plotted.  We  show  here  only  the 
curves  for  the  4%  -  96%  Ar  mixture.  Figure  4a  shows  the  temperature 

behind  the  shock,  and  Fig.  4b  gives  the  density  ratio  across  the  shock, 
p^/p i»  60th  plotted  as  functions  of  the  shock  speed,  Us*  The  three  curves 
on  each  graph  give  the  conditions  when  equilibration  of  the  shocked  gas 
has  proceeded  to  the  states:  i)  translational  and  rotational  equilibrium 
only,  ii)  translational,  rotational,  electronic  and  vibrational  equilibrium 
only,  and  iii)  full  equilibrium  -  but  without  ionization.  In  Fig.  4b,  the 
increase  of  P^/Py  a*  large  Ug  for  the  "Full  Eq"  curve  is  due  to  electronic 
excitation  of  the  atoms.  The  "Vib  Eq"  curve  actually  crosses  the  "Full 
Eq"  curve  due  to  the  fact  that  more  energy  is  invested  in  electronic  exci¬ 
tation  of  than  in  20.  The  three  sets  of  curves  should  not  be  construed 
as  showing  distinct  stages  in  the  relaxation  processes  in  the  shocked  gas 
but  they  do  give  an  indication  of  the  range  of  variation  of  T  and  p  during 
the  relaxation. 


For  the  conditions  existing  in  the  shock -heated  gas  in  these 

experiments,  the  argon-electron  collision  frequency  is  much  less  than 

the  microwave  frequency,  and,  the  shock-heated  plasma  behaves  like  a 

pure  dielectric.  Under  these  conditions  the  electron  density  necessary  to 

#  4 

produce  complete  reflection,  n  ,  is  given  by 


7Tm 


n 


f‘ 


(5) 


where  m  is  the  electron  mass,  f  is  the  microwave  frequency  and  e  is 

the  electron  charge.  For  f  =  1.  5  x  10  sec"  ,  ng  =  2.  8  x  lO1^  electrons/cc 

at  the  lowest  densities  employed  in  this  work,  this  corresponds  to 

n/n.  =8x10^,  which  is  thermally  insignificant, 
c  Ar 

The  microwave  reflection  probe  was  calibrated  with  copper 
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sheets  to  find  its  response  to  a  perfect  reflector  moving  across  the  width 

of  the  quartz  window  in  the  knife-edged  plate  at  the  end  of  the  wave  guide.  In 

2 

Fig.  5  the  reflection  coefficient,  /R/  ,  defined  by  the  equation 

/  /2  _  power  reflected  . 

'  '  power  input  ’ 

is  plotted  vs.  the  fraction  of  the  quartz  window  covered  by  the  copper 

sheet.  Two  different  thicknesses  of  copper  were  used  in  this  calibration, 

both  yielding  essentially  the  same  results.  Note  that  when  the  quartz 

2 

window  is  1/2  covered  by  the  copper  reflector,  /R/  =  .  85. 

When  the  microwave  apparatus  is  mounted  bn  the  shock  tube, 

the  quartz  window  is  bisected  by  the  plane  defined  by  the  slits  in  the  U.  V. 

apparatus.  Consequently,  when  the  U.  V.  signal  first  shows  the  arrival 

of  the  shock  wave,  the  microwave  window  is  already  1/2  covered  by  the 

shocked  gas.  Defining  t  as  the  time  (subsequent  to  the  shock  front 

e  2 

arrival  at  the  U.  V.  slits)  when  the  microwave  signal  /R/  =  .  85,  i.  e.  , 

when  ne  ^  ng  ,  it  is  clear  that  if  TdAe  <  1  no  significant  ionization  has 
occurred  within  the  chemical  relaxation  zone.  Inspection  of  the  dual  trace 
oscillograms  which  display  both  the  U.  V.  and  microwave  signals  showed 
that  this  condition  is  satisfied  in  all  cases.  (For  examples,  see  the 
oscillograms  shown  in  Figs.  6,  7  and  8.) 

The  transmission  of  radiation  of  intensity  Iq  by  the  shock- 
heated  oxygen  is  given  by  the  following  equation: 

I/IQ  =  exp  j^-ki  (p2/Pl)  (P/760)  (02)  j  ,  (7) 

where  I/Iq  is  the  fraction  of  incident  radiation  transmitted,  k  is  the 
absorption  coefficient  per  cm  of  path  length  ,  p2/p^  *8  density 
ratio  across  the  shock  front,  is  the  initial  pressure  in  the  shock  tube 
(mm  of  Hg),  and  (O/  is  the  fraction  of  oxygen  in  the  gas  mixture  at  any 
time.  For  the  mixtures  studied  in  these  experiments  the  temperature 
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of  the  shock-heated  gas,  T2,  and  p2/pj  are  essentlally  constant  behind 
the  shock.  Hence,  once  vibrational  equilibrium  is  obtained  k  is  also  a 
constant.  We  can  therefore  write 


(Oz)  =  -  (1/Q)  In  I/I0 


(8) 


where 

Q=k^£(p2/pi)  (Pj/760)  =  constant.  (9) 


We  can  write  the  equation  for  the  rate  of  disappearance  of  C>2  as 

follows 


d  In  (02)  d  In  (02) 

"  (p2/P1)dtL  “  "kd  (nAr/N0) 


(10) 


where  t  and  tT  are  particle  and  laboratory  times,  respectively.  Here 
we  have  neglected  the  contribution  made  by  C?2  -  O 2  and  O 2  -  O  collisions 
due  to  their  low  frequency  of  occurrence  (see  Sec.  IV  for  a  discussion  of 
this).  Also,  we  have  omitted  the  recombination  term  since  we  do  not 
follow  the  process  close  enough  to  equilibrium  to  make  this  term  signi¬ 
ficant. 

Combination  of  Eqs.  (8)  and  (10)  leads  to  the  following: 


d  log  (log  IQ/I) 


dt. 


~d 

m 


^nAr/N0^ 


(ID 


A  semi-log  plot  of  log  I_/l  vs,  t  should  be  a  straight  line  whose 

u  Xj 

slope  yields  k^.  Figures  6-9  show  four  such  plots  along  with  the  oscil¬ 
lograms  from  which  they  were  produced.  As  is  seen,  the  data  yield 
straight  lines  when  plotted  in  this  manner. 

It  is  to  be  noted  that  the  initial  signal,  L,  (before  shock  wave 
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0.01 
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Fig,  6  Log  Iq/I  plotted  vs.  laboratory  time  for  a  particular  run.  The  inset 
shows  the  oscillogram  from  which  this  plot  was  prepared.  Indicated 
in  the  figure  are  the  resolution,  t  ,  and  the  vibrational  relaxa¬ 
tion  time,  rv,  obtained  by  extrapo5aftng  the  results  of  Ref.  2.  For 
other  notations,  see  text.  Shock  conditions:  4%  0?,  P,  =  .263  mm, 
Us  =  4.  41  mm/jisec,  T  =  16440°K. 
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Fig.  7  Shock  conditions:  4%  O,,  P.  =  .236  mm,  U  =3.53  mm /u 
T  =  10545°K.  2  1  8 
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Fig.  8  Shock  conditions:  2%  O,,  P,  =  1.054  mm,  U  =3.05  mm/(isec, 
T  =  8220°K.  L  1  8 
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Fig.  9  Shock  conditions:  1/2%  O,,  P,  =  10.09  mm,  U  =2.56  mm/usec, 
T  =  6075°K.  L  1  8 
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arrival)  has  already  been  attenuated  by  the  amount  F/Iq.  For  all  the  runs 
above  7700°  (4%  C>2  and  2%  C^),  the  relaxation  times  were  fast  enough 
so  that  a  measurement  of  Iq  could  be  made  subsequent  to  dissociation. 

For  these  runs  this  Iq  was  used  in  reducing  the  data.  When  the  experi¬ 
mental  value  of  F/Iq  was  compared  to  the  theoretical  value  calculated 
from  Eq.  (7)  (setting  k  =  315  and  P^/p^  =  l)  the  agreement  was  always  within 
a  few  percent.  For  the  low  temperature  runs  (1/2%  O^),  the  relaxation 
times  were  too  long  to  follow  the  signal  to  Iq,  i.  e.  ,  boundary  layer  build¬ 
up  and  change  of  light  source  intensity  made  this  too  uncertain.  In  these 
cases,  the  theoretical  F/Iq  was  utilized' to  obtain  Iq. 

Indicated  in  Figs.  6-9  are  resolution  times,  r  ,  which  were 
computed  by  the  equation 


T 

res 


S  +  6  /4 
- D - 


(12) 


where  the  slit  width  S  =  1  mm,  and  6  is  the  shock  front  curvature  given 
by  Eq.  (4).  This  theoretical  resolution  agrees  well  with  the  experi¬ 
mentally  determined  signal  rise  time,  the  exception  to  this  being  the 
=  4  and  10  mm  runs  (l/2%  C^)  where,  because  of  the  high  density, 
there  was  an  optical  effect  at  early  times  probably  associated  with  internal 

30C 

reflections  off  the  shock  front.  This  latter  effect  persisted  in  pure  Ar 

shocks  run  at  comparable  densities.  For  this  reason,  these  oscillograms 

were  not  read  at  times  less  than  2  /j,sec. 

Also  indicated  on  these  figures  are  the  vibrational  relaxation 

times  T  computed  from  Eq.  (3)  for  the  experimental  conditions.  The 

condition  t  /t  >  1  is  met  in  these  experiments  only  at  low  densities 
v'  res 

and  temperatures,  in  particular,  for  the  4%  runs  at  the  lower  tempera- 


5  A  similar  effect  to  this  was  observed  in  previous  work  at  12 70 A, 
see  Ref.  5. 
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tures  at  which  this  mixture  was  studied.  For  example,  see  Fig.  7  -  notice 
how  the  data  deviate  from  the  straight  line  at  t  <  r^.  Under  these  condi¬ 
tions,  the  absorption  coefficient,  k,  can  no  longer  be  treated  as  a  constant 
in  Eq.  (7),  consequently,  Eq.  (11)  is  no  longer  valid.  The  upper  point 
located  on  the  zero  time  axis  of  Figs.  6-9  and  labeled  "No  Vib"  was  cal¬ 
culated  from  Eq.  (7)  setting  k  =  315  and  using  the  density  ratio  appropriate 
to  no  vibration.  Notice  how  the  data  tend  toward  this  origin  in  Fig.  7  for 

t  <  t  . 

L  v 

The  theoretical  t^  =  0  intercept  for  these  straight  line  plots  can 
be  computed  via  Eq.  (7)  using  the  value  of  k  chosen  consistent  with  complete 
vibrational  equilibrium  and  p ^/P\  chosen  for  the  state  "vibration  only". 

This  was  done  for  all  the  runs.  For  Figs.  6-9  these  points  are  identified 
as  "Vib".  In  all  except  some  of  the  lowest  temperature  runs  (below  6500^) 
the  extrapolated  straight  line  given  by  the  data  intercepts  the  ordinate  above 
the  theoretical  point.  This  implies  that  the  dissociation  of  does  not 
start  at  zero  time  with  the  rate  that  is  ultimately  measured.  Indeed,  there 
is  an  incubation  time,  At,  during  which  essentially  no  dissociation  occurs, 
this  time  being  given  by  the  intersection  of  the  extrapolated  straight  line 
with  the  "Vib"  ordinate. 

The  zero  of  time  was  chosen  from  the  oscillograms  at  the 
beginning  of  the  signal.  The  resolution  in  this  experiment  was  limited  by 
the  combination  of  finite  slit  width  and  shock  front  curvature.  Because  the 
experimental  measurement  of  I/Iq  averaged  over  this  finite  slit  and  contin¬ 
uously  looked  through  a  gas  sample  which  had  a  radial  time  displacement 
due  to  the  shock  front  curvature,  a  slightly  displaced  time  scale  results. 
Using  the  rate  constant  deduced  from  each  run,  the  shock  conditions,  the 
slit  width  and  the  shock  front  curvature,  a  theoretical  curve  waq  computed 
on  an  IBM  650  computer  which  integrated  out  the  effect  of  the  finite  slit 
and  shock  front  curvature.  The  pertinent  equations  for  this  procedure  are 
derived  in  the  Appendix.  The  experimental  curve  would  be  identical  with 
this  computed  curve  if  the  experimental  curve  actually  started  off  with  its 
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ultimate  slope.  These  curves  were  plotted  alongside  the  experimental 
curves;  they  were  straight  lines  having  the  same  slope  as  the  experimental 
curve  but  the  intercept  of  this  computed  curve  fell  between  the  intercept 
of  the  experimental  curve  and  the  theoretically  computed  intercept  which 
did  not  take  into  account  shock  front  curvature  and  finite  slit  width.  The 
time,  At,  at  which  the  extrapolated  experimental  line  crossed  the  intercept 
ordinate  of  the  computed  curve  represents  a  time  during  which  the  dissocia¬ 
tion  process  is  inhibited.  These  intercept  ordinates  are  labeled  "Vib  (Cor¬ 
rected)"  in  Figs.  6-9. 


IV.  Results  and  Discussion 

An  Arrhenius  plot  of  the  dissociation  rate  covering  the  tempera¬ 
ture  range  5000  -  18000°K  was  constructed  and  is  shown  in  Fig.  10.  The 
temperatures,  T,  in  this  figure  are  the  average  of  the  "Vib  Eq"  and  the 
"Full  Eq"  temperatures.  The  least  squares  fit  to  the  data  between  5000 
and  11000°  yielded  an  -  Ar  dissociation  rate  given  by  the  following 
equation 

kd  =2.9  (+  12%)  x  1014  exp  (-D/RT)  cc/mole-sec,  (13) 

in  which  the  dissociation  energy  of  is  D  =  118000  cal/mole.  The  fit  to 
this  equation  was  carried  out  by  computing  the  value  of  k^/exp  (-D/RT) 
for  each  data  point  and  then  taking  the  average  for  all  the  data  points. 

The  probable  error  was  computed  from  the  standard  deviation.  Although 
the  data  could  have  been  fit  with  an  equation  having  a  pre -exponential 
temperature  factor  in  it,  as  other  investigators  have  done,  it  was  not  felt 
that  the  data  warranted  such  treatment. 

At  temperatures  above  11000°  the  data  points  fell  below  the 
Eq.  (13)  extrapolated  straight  line,  the  deviation  growing  larger  at  higher 
temperatures.  At  18000°  k^  is  about  0.  45  times  the  value  given  by  the 
above  equation.  Also  shown  in  Fig.  10  are  the  results  obtained  by  Camac 
and  Vaughan  [Eq.  (1)]  which  are  extrapolated  for  T  >  7500°.  In  a  recent 
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Fig.  10  The  rate  constant  for  the  dissociation  of  O,  by  Ar  plotted  vs. 

reciprocal  average  temperature.  The  dasned  line  is  from  the  work 
of  Camac  and  Vaughan  and  is  extrapolated  for  T  >  7500°K.  The 
least  squares  fit  to  the  data  shown  by  the  solid  line  yields  k  ,  =  2.  9 
(+  12%)  x  10*4  exp  (-D/RT)  cc/mole-sec.  d 
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report,  Anderson,  ^  in  an  investigation  similar  to  that  of  Camac  and 
Vaughan,  gives  the  02  -  Ar  dissociation  rate  for  temperatures  between 
4000  and  6000°  as 

IO  Q  Q 

kd  =  l.  80  xl0AO  T  *7  exp  (-D/RT)  cc/mole-sec.  (14) 

This  yields  a  rate  constant  1/2  as  large  as  that  of  Camac  and  Vaughan  and 
therefore  is  in  excellent  agreement  with  the  present  data  where  the  tempera¬ 
tures  overlap  (5000  -  6000°). 

As  was  already  mentioned,  02  -  02  and  C>2  -  O  collisions  were 
neglected  in  this  work.  The  results  of  several  independent  investigations* 
indicate  that,  at  least  at  low  temperatures 

kd  (Oz  -  Ar)  =  i  kd  (02  -  02)  =  i  kd  (02  -  O)  .  (15) 

The  low  temperature  work  (T  <  7700°)  in  this  investigation  was 
carried  out  with  the  1/2%  C>2  mixture.  Concentrating  our  attention  on  the 
C>2  -  O  collision  (which  would,  of  course,  have  the  largest  effect)  we  find 
that  the  small  increase  in  the  effective  rate  constant  caused  by  the  forma¬ 
tion  of  O-atoms  (a  factor  of  «  1.  25)  is  almost  exactly  cancelled  by  the 
decrease  in  the  rate  constant  caused  by  the  corresponding  temperature 
drop  («115°).  On  the  other  hand,  for  the  2%  and  4%  mixtures  at  high  tem¬ 
peratures  one  should  observe  an  effective  rate  constant  which  increases 
with  time  if  the  relations  of  Eq.  (15)  apply.  For  example,  for  the  case 
4%02  and  T  =  17000°K,  the  O-atoms  should  increase  the  effective  rate 
constant  by  a  factor  of  «3  by  the  time  the  02  is  all  dissociated;  the 
decrease  in  the  rate  constant  due  to  the  750°  temperature  drop  only  yields 
an  opposing  factor  of  0.9.  Since  the  data  yielded  straight  log  Iq/I  vs.  t^ 


?  These  are  discussed  and  summarized  in  Ref.  7. 
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plots  for  the  4%  o 2  cases  at  high  temperature,  it  is  clear  that  the  effective 
rate  constant  is  essentially  constant  during  the  dissociation  process. 
Therefore,  we  conclude  that  in  the  vicinity  of  17000°  kd  (C>2  -  O)  »1.  3  kd 
(Oz  -  Ar).  This  is  not  too  surprising,  for  at  17000°  the  product  of  the 
collision  number  and  the  chemical  relaxation  time  is  234,  i.e. ,  the  average 
oxygen  molecule  dissociates  in  234  collisions  with  Ar.  To  maintain  the 
catalytic  efficiency  ratio  of  25  at  17000°  implies  that  it  would  take  but 
9  Oz  -  O  collisions  to  dissociate  the  average  O^.  It  is  clear  that  this 
factor  of  25  must  decrease  at  high  temperatures. 

The  uncorrected  incubation  times,  all  reduced  to  a  common 
17 

density  (1.28  x  10  particles/cc),  are  plotted  vs.  Ty.^  in  Fig.  11.  It 
is  to  be  noted  that  below  6500°  some  measured  incubation  times  were 
negative  and  are  plotted  in  this  figure  as  solid  symbols.  One  point  at 

T  =  6500,  At_  =  0  is  not  plotted. 

Vib  Eq  P 

In  Fig.  12  the  same  data  are  plotted  but  now  corrected  for  the 
effect  of  the  finite  slit  width  and  shock  front  curvature  as  outlined  in  the 
Appendix.  We  have  omitted  the  runs  below  6500°  for  =  4  mm  and  below 
5500°  for  Pj^  =  10  mm  because  of  the  tremendous  scatter  in  that  data.  (See 
discussion  of  errors  below.  )  Comparison  of  Figs.  11  and  12  shows  that  the 
corrections  to  At  (decreasing  the  apparent  incubation  time)  were  largest 
at  the  higher  temperatures  as  one  would  expect.  As  seen  in  Fig.  12,  the 
incubation  time  decreases  by  about  two  orders  of  magnitude  between  6500 
and  18000°. 

In  Fig.  13  we  have  plotted  the  ratio  of  the  corrected  incubation 
time  to  the  vibrational  relaxation  time  obtained  from  Eq.  (3).  The  ratio 
At/r  is,  of  course,  density  independent.  The  solid  line  drawn  through 
the  data  in  Fig.  13  is  somewhat  subjective;  however,  it  will  serve  to  allow 
some  comparisons  to  be  made. 


Alp  “UNCORRECTED"  ( fiSEC ) 


Fig.  11  The  uncorrected  incubation  times  in  the  particle  coordinate  system 
at  a  density  of  n^r  =  1.28  x  101'  particles/cc  plotted  vs.  the  "Vib 
Eq"  temperature. 
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Let  us  define  the  degree  of  vibrational  excitation,  r)  ,  necessary 
for  an  uninhibited  dissociation  rate  by  the  ratio 


Tl  (Tt)  = 


(16) 


where  Tv  and  Tt  are  the  vibrational  and  translational  temperatures  respec¬ 
tively.  ti  is  itself  a  function  of  translational  temperature. 


The  vibrational  energy  relaxes  according  to  the  equation 

E 


8 


v 

T 


=  exp  ( -t/T  v) 


(17) 


where  E^  is  the  vibrational  energy  corresponding  to  the  final  conditions 
and  Ev  is  the  vibrational  energy  at  any  time,  t. 

Ev/Ej  is  given  by  statistical  mechanics  as 

E_  exp  (6/Tf)  -1 

expTS/T* )  -T 


v 

"*7 


(18) 


where  is  the  final  vibrational  temperature  consistent  with  E^  and  Tv  is 


the  vibrational  temperature  corresponding  to  E 


We  wish  to  calculate  T  /T.  = 

v'  t 


For  02,  6  =  2228  K. 


when  the  time  t  =  At.  Noting 
is  identical  with  the  translation  - 


that  the  final  vibrational  temperature,  T^, 
al  temperature  Tt>  we  combine  Eqs.  (16),  (17)  and  (18)  to  give 

At  ,  (  exp  (  0 /t\  Tt)  -1  ) 

—  n  \  exp(07TlTt)  -«p^/Tt)  ] 


(19) 


In  Table  I  we  give  T)  calculated  from  Eq.  (19)  at  several  tempera¬ 
tures.  The  At/r  at  each  temperature  corresponds  to  the  solid  line  in 
Fig.  13. 
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Table  I 


Tt 

A  t/r 

V 

3 

11 

5500 

2.4 

.93 

8000 

1.  1 

.71 

11000 

.  55 

.  48 

14000 

.  45 

.  41 

18000 

.40 

.  37 

/T 


t 


Although  the  scatter  in  the  data  is  considerable,  the  results  clearly  indicate 
that  the  dissociation  rate  of  C>2  is  severely  limited  until  some  degree  of 
vibrational  equilibration  obtains,  the  degree  of  equilibration  necessary 
being  higher  at  lower  translational  temperatures. 

The  largest  source  of  error  in  this  work  is  due  to  the  uncertainty 
in  IQ.  This  is  caused  by  a  combination  of  reading  errors  on  the  oscillo¬ 
grams,  boundary  layer  effects,  and  light  intensity  variation. 

The  straight  line  of  Eq.  (11)  can  be  represented  by 

y  =  log  (log  Iq/I)  =  a  t  +  b  .  (20) 


Eq.  (20)  yields 


At  =  (  yQ  -  b  )/a  .  (21) 


where  yQ  is  the  ordinate  at  t  =  At.  For  the  two  coordinates  (y2>  t2)  and 

(yL*  tj)  we  have 


a  = 


y2  -n 

t2  "  h 


and 


b  = 


yi  t2  -  y2  tl 

h  '  h 


(22) 


Differentiation  of  the  above  equations  yields  the  following  relations; 


Ay  =  (.189)  10"y 


(23) 
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Aa 

a 


A  y2  -  AyL 

'fz  '  n 


Ab_  _  AVj  t2  -  Ay2  tl 

a  ^2  -n 


A  (At)  _  A  b  A  a 

51  -  "a  At  '  a 


(23) 

(Cont.  ) 


With  the  aid  of  these  equations,  the  fractional  errors  A(At)/At  and  Aa/a 
were  computed  at  several  temperatures  using  typical  values  of  At,  y^,  y2> 
t.  and  t_  consistent  with  the  experiment.  These  errors  are  tabulated  in 
Table  U.  We  have  taken  A Iq/Iq  to  be  4.  0  x  10  ~  ,  which  is  probably  an 
upper  limit  to  the  uncertainty  in  A  Iq/Iq ‘ 


Table  II 


A  Iq/Iq  assumed 

+  0.  04 

T°K 

P^(mm) 

A(At)/At 

Aa/a 

5500 

10 

+  L  3 

T  .052 

6000 

10 

+  .  58 

T  .094 

6000 

4 

+  2.2 

T  -14 

8000 

1 

+  1.  0 

+  -55 

12000 

1/2 

+  .  53 

T  *68 

17000 

1/4 

±  *14 

00 

'O 

1+ 

These  error  calculations  are  the  bases  for  discarding  the  At 
data  below  6500°  for  the  4  mm  runs  and  below  5500°  for  the  10  mm  runs. 
Notice  that  these  error  calculations  indicate  a  larger  fractional  error 
in  the  rate  constant  (Aa/a)  at  the  higher  temperatures.  Inspection  of  Fig.  10 
bears  this  out. 
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APPENDIX 


Here  we  derive  the  equation  analogous  to  Eq.  (11)  but  embodying 
corrections  for  the  shock  front  curvature  and  finite  slit  width.  The  follow¬ 
ing  sketch  will  be  necessary  to  define  the  problem. 


The  shock  front  (curvature  exaggerated  in  this  sketch)  has 
moved  across  the  finite  slits  of  width  S  (exaggerated  width)  from  left  to 
right.  Time  equals  zero  when  the  shock  arrived  at  the  position  x  =  y  =  0. 
Two  gas  elements  are  shown  by  the  dashed  curves.  It  is  clear  that- 

(R  -  x')2  +  y2  =  R2  ,  (A-l) 

2 

and,  ignoring  the  small  (x1)  term, 

x'  =  y2/2R  ;  6  =  (_^/2)2/2R  ,  (A-2) 

where  the  6  is  understood  to  be  the  shock  front  curvature  existing  between 
the  knife  edges. 


We  define  the  time  t'  to  be  the  time  elapsed  between  the  arrival 
of  the  shock  front  (along  y  =  0)  at  the  infinitely  narrow  slit  and  the  arrival 
of  the  gas  element  #2;  i.  e. , 

f  =  t  -  x/Ug  •  (A-3) 

It  is  clear,  however,  that  the  time  the  gas  has  been  heated  decreases  as  we 
move  out  from  y  =  0  to  y  =  ^ /l\  indeed,  at  the  point  (x,  y)  we  are  looking 
at  gas  (in  gas  element  #1)  which  has  been  heated  for  a  time  given  by 


t"  =  t'  -  x'/U 

8 


(A-4) 


This  is  perfectly  general  and  t"  is  the  time  which  appears  in  the  exponen¬ 
tial  decay  law  for  the  dissociation.  We  can  now  write 

exp  (-t'/r) 

where  t  is  the  dissociation  relaxation  time  defined  by  Eq.  (2)  but  in  the 
laboratory  coordinate  system. 

The  absorption  equation  [Eq.  (7)]  written  in  differential  form  is 


exp  (y  /2R  U  r  )  ,(A-5) 


']■' 


(07)  =  exp 


[- 


(f  -  x'/Us)/r 


] 


d  In  I/I0  =  -  Q  (02)  dy 


(A-6) 


where  Q  is  defined  in  Eq.  (9).  Elimination  of  (02)  between  (A-5)  and  (A-6) 
and  integration  across  the  optical  path  from  / 2  to  +^/2  yield 


=  exp 


2  Q  exp  ( -t'/ r  ) 


i 


J!/  2 


exp  (y  /2  R  U8  T  )  dy 


(A-7) 
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Making  the  substitution 


y*  =  y/<2R  Ufl  r)  l/Z 
and  with  the  aid  of  (A-2),  we  find  (A-7)  becomes 

l/l0  =  exp  £  -Q  Jl  p  exp  (  -t'/T ) 
\l/2 


where 


-  ■  m'"/ 


■w) 


exp  (y1)  dy* 


The  next  step  is  to  average  over  the  finite  slit,  hence  we 


/  i  t  ■  / 


exp 


-  q/p  exp(  -t '/  t  ) 


dx 

S 


and  making  the  substitution 


Z  =  x/U  r 
'  8 


we  have 


t  -ft-')/ 


s 

U~T~ 


exp 


-  P  exp(  -t/ t  )  exp(Z) 


dZ 


Defining 


y  =  P  exp  (-t/T) 


(A-8) 


(A-9) 


(A-10) 


write 

,  (A-ll) 


(A -12) 


.  (A-13) 


(A-14) 
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which  is  independent  of  x  ,  and  making  the  substitution 


U  =  y  exp  Z 


we  find 


exp  (S/Us  t  ) 
exp  (  -U  ) 


dU 

U 


(A-15) 


(A-16) 


For  a  given  shock  tube  run,  6/U  t  is  known;  hence  p  can  be 
evaluated  from  (A-10)  and  is  independent  of  time,  and  -y  can  be  evaluated 
from  (A-14)  as  a  function  of  time  and  the  corrected  i/l ^  can  be  calculated 
from  (A-16)  as  a  function  of  time. 
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